Probing weakly-bound molecules with nonresonant light 
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We show that weakly-bound molecules can be probed by "shaking" in a pulsed nonresonant laser 
field. The field introduces a centrifugal term which expels the highest vibrational level from the 
potential that binds it. Our numerical simulations applied to the Rb2 and KRb Feshbach molecules 
indicate that shaking by feasible laser pulses can be used to accurately recover the square of the 
vibrational wavefunction and, by inversion, also the long-range part of the molecular potential. 

PACS numbers: 37.10.Vz, 34.20.Cf, 34.50.Ez, 33.90.+h 



Weakly-bound species, such as haloes [lj], spend most 
of their vibrational period in the classically forbidden re- 
gion of the electronic potential. The accuracy required 
to treat such a potential stretches current computational 
methods to their limits. Here we present a versatile tech- 
nique for determining the interatomic potential from the 
dependence of the dissociation probability of the weakly- 
bound molecule on the intensity of a pulsed nonresonant 
laser field. We also show that a cw laser field can be used 
to control the molecules' size. 

Since diatomic molecules formed by photo- or 
Feshbach-association of ultracold atoms often emerge 
in the highest vibrational level of a ground electronic 
state @-|U, they represent prototypical quantum halo 
systems. Among other prominent examples of haloes arc 
the recently observed Efimov trimer [B[ and tetramer 6 1 , 
the latter produced in collisions involving halo dimers 7] . 

The technique relies on the ability to impart a value of 
angular momentum to the weakly-bound molecule such 
that the centrifugal term concomitant with it expels the 
molecule's vibrational level from the potential and thus 
causes the molecule to dissociate. Unlike in the case of 
rotational predissociation, the value of the imparted an- 
gular momentum is tunable, derived from the anisotropic 
polarizability interaction of the molecule with a nonres- 
onant laser field. This interaction creates directional 
states in which the molecular axis librates (shakes) about 
the field vector The laser intensity needed to im- 
part a preordained value of the angular momentum varies 
characteristically with the internuclear distance. It is 
this characteristic dependence that can be used to map 
out the probability density of the vibrational state from 
which the molecule was forced to dissociate. A highly 
accurate long-range molecular potential can be then ob- 
tained by inverting the vibrational probability density. 
This route to an accurate potential, independent of spec- 
troscopy or scattering, complements what can be learned 
from either. We illustrate the technique's machinery by 
examining Feshbach molecules of acute interest, Rb2 and 
KRb. 

Although the technique is applicable to weakly-bound 
molecules in any rotational state, we limit our con- 
siderations here to the case of rotationless (diatomic) 



molecules, i.e., species in vibrational levels that can only 
support the ground rotational state. In addition, we con- 
sider the molecule to be in a 1 E electronic state, with a 
V(r — > oo) = — Cer~ e asymptotic potential. In order for 
a level to be rotationless, its binding energy, Eb, must 
fulfill the inequality \Ef,\ < d^m" 3 / 2 ^ 1 ^ 2 , where m is 
the reduced mass and d§ w 1.6 is a dimcnsionless param- 
eter which depends solely on the exponent of the power 
law potential and can be evaluated analytically Q. The 
centrifugal term in the effective potential 



U(r) = V(r) + 



{j 2 )h 2 

2mr 2 



(1) 



then leads to dissociation for any higher rotational 
states. Here V(r) is the electronic potential and (J 2 ) = 
(tp{t)\ J 2 IV'(^)) is the expectation value of the square of 
the imparted angular momentum in a state \iMt)) created 
by the anisotropic polarizability interaction [9| . The crit- 
ical value of the angular momentum, J*, corresponding 
to (J* 2 ) = J* (J* + 1), for which the effective poten- 
tial can no longer support a given vibrational level can 
be much smaller than unity. For instance, for the high- 
est vibrational state, v = 123, of 85 Rb 2 ( 1 £+), we find 
J* = 0.22. Such an angular momentum can be imparted 
to the molecule by a feasible laser field. As noted earlier 
[icj |. for 4 He2( 1 £g ,v — 0), the critical angular momen- 
tum J* is 0.0337. 

We obtain the wavefunction \ip(t)) by solving the 
time-dependent Schrodinger equation for a dimensionless 
Hamiltonian 



Hit) 
B 



g(t) [Aw cos 2 9 - 



(2) 



Here g(t) is the time profile of a plane-polarized laser 
pulse of peak intensity /, Auj = u>\\ — lu±, and j_ = 
2iva\\ t ±I/(Bc), with Q!||(r) and a±(r) the polarizability 
components parallel and perpendicular to the molecular 
axis, and B = B(r) = h 2 /(2mr 2 ). Because of the az- 
imuthal symmetry about the field vector, the induced 
dipole potential involves just the polar angle 9 between 
the molecular axis and the polarization plane of the laser 
pulse. We assume the oscillation frequency of the laser 
field to be far removed from any molecular resonance 
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FIG. 1: Temporal dependence of the angular momentum im- 
parted to a molecule by laser pulses of duration r = 2.5 irh/B 
(green solid line), r = 0.25 nh/B (blue solid line), and 
r = 0.025 irh/B (red solid line). The corresponding time pro- 
files of the laser pulses, centered at time 3nh/B, are shown 
by dotted lines. Also indicated is the angular momentum due 
to a cw laser field (black solid line). Note that Auj = 100 at 
the peak of the pulsed fields and for the cw field. See text. 



and much higher than the reciprocal of the pulse du- 
ration; as a result, the time dependence of the radia- 
tive field is reduced to that of the time profile [ijj . For 
homonuclear molecules, the dependence of the polariz- 
ability anisotropy, Aa(r) = cei i (r) — a±(r), on the in- 
ternuclear distance is well captured at large r by Sil- 
berstein's expansion, Aa(r) = &a^r~ 3 + 6agr -6 + 
where ao is the atomic polarizability |12j. The solutions 
of the reduced time-dependent Schrodinger equation cor- 
responding to Hamiltonian @, 



. h &ip(t) _ H(t) 
B dt ~ B 



(3) 



can be expanded in a series of field-free rotor wave func- 
tions \JM) = Yjm (pertaining to eigenenergies Ej) 



^(Aw(t)) = ^ej(Au;(i))|J,M) ex P 



iEjt 



(4) 



whose time-dependent coefficients, cj (Aw (t)), solely de- 
termine the solutions at given initial conditions (in the 
interaction representation) . Thus the hybridization coef- 
ficients cj can be found from the differential equations 



h 



i— cj (t) = ~y~]cj>(t) (J,M\ Au> cos 2 6 + uj±_ \ J',M) 
j' 

i(Ej,-Ej)t 



x exp 



g(t) (5) 



where we consider the pulse shape function to be a Gaus- 
sian, g(t) = exp [— 4 ln(2)t 2 /r 2 ] , characterized by a full 
width at half maximum, r, the "pulse duration." By 
taking into account the non-vanishing matrix elements 
(J, M\ cos 2 9 | J', M), Eq. (0 reduces to a tridiagonal 
form and can be numerically solved, at any particular 



value of r, by standard methods We note that 

the expectation value of B in a given vibrational state, 
(j) v = \v), is the rotational constant pertaining to that 
state, B v = (v\B\v) = -^(v\r~ 2 \v) . Similarly, the ex- 
pectation value Aa v = {v\Aa\v} is the polarizability 
anisotropy of the molecule in the vibrational state \v). 

Figure[T]shows the temporal dependence of the angular 
momentum imparted to a molecule at a time for 
different values of the pulse duration r. One can see that 
whereas for r > 2.5 irh/B, the angular momentum of the 
molecule is enhanced only during the time the laser pulse 
is on (adiabatic limit), for r < 0.25-7T h/B the molecule 
acquires a part of the transferred angular momentum for 
good (nonadiabatic limit). At r < 0.025 nh/B, most of 
the angular momentum is imparted for good. Note that 
the amplitude of the laser pulse, Aoj — 100, is the same 
for all three pulses. The same value of Aui was also used 
for the cw, stationary case, with g(t) = 1, shown by the 
horizontal line. 

Thus, in the nonadiabatic regime, characterized by 
r <g; nh/B, a laser pulse can endow a molecule, on a 
timescale r, with a preordained value of angular mo- 
mentum which the molecule will keep after the pulse 
has passed. For pulses shorter than the vibrational pe- 
riod, t < t v , the interaction with the laser field will be 
well within the nonadiabatic regime. Moreover, in this 
case, the imparted angular momentum will also depend 
on the value of the internuclear distance, r — r pi which 
the molecule had at the moment when the laser pulse 
struck. For disparate time scales such that t <^t v <C r r , 
the dimer is transferred from the V(r) to the U(r) po- 
tential essentially instantaneously. 

The laser intensity needed to impart a given value of 
angular momentum at a given internuclear distance can 
be found by solving Eq. |3|). For a critical value (J* 2 ) of 
the angular momentum, Eq. f3| yields a laser intensity 
/ which corresponds to a critical internuclear distance 
r* such that when struck at r p — r* , the molecule will 
dissociate. It is shown below that I(r*) increases essen- 
tially monotonously with r* . As a result, for a given 
intensity I(r*), the molecule will receive a critical value 
of the angular momentum and dissociates if the pulse 
strikes at any distance r p < r* . On the other hand, for 
r p > r* , the imparted angular momentum is subcritical 
and the molecule will remain intact. Hence the dissocia- 
tion probability coincides with the probability of finding 
the molecule with an internuclear distance less than r* , 



\4>v{r)\ dr. 



(6) 



The dependence, F — F(I), of the dissociation proba- 
bility F on the intensity / of the laser pulses is exper- 
imentally observable. By making use of the computed 
/ = I(r*) dependence, it can be transformed into F(r*), 
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FIG. 2: The "exact" single potential energy curve V(r) of 
85 Rb2 (green solid curve) and the wavefunctions of the high- 
est vibrational state, v — 123, obtained from the "exact" 
potential (red solid line) and from the analytic model (blue 
solid line). The binding energy of the v = 123 level is 
E b = -237 kHz. See text. 



whose derivative with respect to r* yields the vibrational 
probability density, |0„(r)| 2 . Note that while F(r*) de- 
pends on the molecular potential, I(r*) does not, as it 
reflects just the molecular moment of inertia and polar- 
izability. In the short-pulse limit, I(r*) has a universal 
asymptotic dependence on r* , I(r* — > oo) cx (r*) 3 , which 
can be derived analytically from the closed-form solution 
of Eq. ^ in that limit, given in ref. (Tj|. The polariz- 
ability anisotropy Aa(r) at large r is accurately captured 
by Silberstein's formula [12|. Alternatively, the need to 
know Aa(r) can be eliminated by measuring F(I) for 
two weakly-bound vibrational levels. 

We now illustrate the outlined procedure by evaluat- 
ing the dissociation probability of the 85 Rb2 molecule 
in its highest vibrational state. Based on the molecu- 
lar constants and a combination of the potential energy 
curve of ref. 



14j with the long-range exchange and dis- 



persion terms of ref. 15[ , a single potential curve for the 
ground 1 S+ electronic state was found, with a dissoci- 
ation energy D e = 3993.53 cm -1 and an equilibrium 
distance r e — 4.21 A. An extrapolation of this "ex- 
act" potential suggests that the highest bound state has 
v = 123. We corroborated this result by numerically 
solving the Schrodinger equation for this potential, which 
also yielded the binding energy and the wavefunction of 
the v = 123 state, see Figure [5] The v = 123 vibra- 
tional state is no longer bound if the squared angular 
momentum is greater than (J* 2 ) = 0.27. We found from 
Eq. ([5]) in the limit t — > oo and for a dependence of 
the molecular polarizability on the internuclear distance 
taken from ref. [1 61 ] that such an angular momentum is 
imparted to the 85 Rb2 molecule by a cw laser field of 
intensity I = 6.58 x 10 8 W/cm 2 . 

In order to map out the wavefunction squared of the 
highest vibrational state, the molecule must be shaken by 
a pulsed laser field whose pulse duration is much shorter 
than the vibrational period. Since the vibrational period 



FIG. 3: Dependence of the intensity of 50 ps laser pulses 
needed for imparting a critical value of angular momentum to 
the Rb2( 1 E + ,u = 123) molecule on the critical internuclear 
distance r*. See text. 



of halo molecules is often in excess of nanoseconds, the 
vibrational motion can be adequately parsed by 50 ps 
pulses, easily available in the laboratory [rjj . The use of 
such pulses simultaneously ensures that the polarizability 
interaction will be well within the nonadiabatic regime, 



as t <c t v <C |p- However, since t v = 0.67 fis for the 
v = 123 state of 85 Rb2, even ns pulses would suffice to 
probe it. 

The calculated dependence of the intensity I of 50 
ps laser pulses needed to dissociate 85 Rb2(t> = 123) 
on the critical internuclear distance r* is presented in 
Fig. [3] The calculation, carried out by solving the time- 
dependent Schrodinger equation ([5]) with (J* 2 ) at each 
distance r = r* , was based on the same molecular pa- 
rameters as the above cw result. Marked in the figure is 
a particular value, /, of the laser intensity which is as- 
sociated with a particular critical internuclear distance, 
f*, as well as the ranges of the internuclear distance r p 
which do or do not lead to dissociation. We note that 
before reaching the asymptotic dependence I oc (r*) 3 
at r* > 10 A, the pulse intensity first decreases and 
passes through a minimum at about r* min — 5.8 A, up to 
which the I(r*) dependence is two- valued. Nevertheless, 
the short-distance region may be included in the calcula- 
tion by simply changing the lower limit of integration in 
eq. ([6]) from zero to some "inner critical distance." How- 
ever, since the time the molecule spends in the r* < r* min 
range is by a factor of about 2 x 10 5 times less than the 
time in spends at r* > r* min , it is a fairly good approxi- 
mation to neglect any effects that may arise at r* < r^ n . 

The vibrational probability density, Fig. [21 combined 
with the dependence of the laser intensity on the critical 
distance I(r*), Fig. [3l yields the dissociation probability, 
Eq. ([6|), as a function of the laser intensity applied, F(I). 
This is shown in Figure [4] by the full curve. Hence a 
measured dissociation curve can be used, by reversing 
the procedure, to reconstruct the vibrational probability 
density and hence the potential that binds it. 

The process of dissociation by shaking is amenable to 
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FIG. 4: The dissociation probability of Rb 2 ( 1 E + ,u = 123) as 
a function of the intensity of 50 ps nonresonant laser pulses. 
Calculation based on the "exact" potential are shown by the 
solid red line, the analytic result by the dashed blue line. See 
text. 



an analytic treatment based on a near-threshold expan- 
sion of the wavefunction [181 ] . The analytic wavefunction, 
obtained for the last bound state of the 85 Rb2 dimer, and 
the corresponding dissociation probability, Eq. ([6|), are 
shown in Figs. [2] and [4] respectively. The oscillations ob- 
served in the exact F(I) are thus identified as due to the 
nodes of the vibrational probability density, and the po- 
sition of the wavefunction's maximum as corresponding 
to the point of maximum slope, "the edge," of F(I). 

We also assessed the dissociation probability of 
39 K 87 Rb in its highest vibrational state, v = 99. Since 
the molecule is polar, we invoked half-cycle pulses, with 
r = 50 ps, to probe it. For a potential taken from 
ref. (l9j and the body-fixed permanent dipole moment de- 
pendence on r from ref. [20j , we found that the molecule 
can be shaken enough to dissociate by laser intensities as 
low as 10 4 -10 8 W/cm 2 for r < 10 A. As the permanent 
dipole moment plummets with increasing r, at higher in- 
ternuclear separations the polarizability interaction pre- 
vails and there is no advantage in using the half-cycle 
pulses. 

We note that a non resonant cw-laser field can be used 
to control the scattering length and positions of magnetic 
Feshbach resonances in a simpler way than recently de- 
scribed by Bauer et al. 



2l|. Furthemore, the laser field 



of optical dipole traps may be able to dissociate some 
of the weakest-bound molecules or to reduce their ap- 
parent binding energy. However, a 100 mW laser beam 
focused to a waist of 50 /im will not, as it imparts only 



a negligible angular momentum of (J 2 ) 



io- 



to a 



S5 



Rb2 molecule. On the other hand, a cw field may 
be used to control the "size" of a halo molecule, i.e., 
(r v ) = (<j) v (r)\r\(j) v (r)), by uplifting the highest vibra- 
tional state v. For instance, 85 Rb2(w = 123), which has 
a field- free mean radius (r v ) — 167 A, is expected to be- 
come a near-threshold state at an intensity of 6.52 x 10 8 
W/cm 2 (corresponding to J = 0.217), with a mean ra- 
dius of 1320 A. This represents an alternative to control- 



ling molecular size by changing the magnetic field in the 
vicinity of a Feshbach resonance Q ■ Conversely, the non- 
resonant field can be used to provide a tunable barrier to 
molecule formation. 
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